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SUMMARY 


A study was made of the thermal radiation associated with the ablation of 
plastic materials . As ablation product materials diffuse into the boundary 
layer, they become possible important sources of radiant energy. This radiation 
may be due to either the increase in temperature of the ablation products or to 
chemical reactions occurring in the boundary layer. 

For the tests, small blunt-nosed projectiles were launched at high veloc- 
ities past photoelectric devices which recorded the resulting radiation. A 
simple mathematical model was developed to describe the boundary injection- 
diffusion process, and equations were derived for the surface heating from this 
source for vehicles of spherical geometry. The test results were used to 
evaluate empirical constants in the equations. 

Three ablation materials were tested: lexan, a polycarbonate resin; General 

Electric Company Century Series 124, a char -forming plastic; and a high-density 
polyethylene . 

Predictions of the stagnation-point heating of large-sized vehicles based on 
the above tests and mathematical model indicate that the radiative heating of 
ablation products varies over a wide range depending upon the ablation material. 
For at least one material the radiant heating is of the same magnitude as the 
convective heating. 


INTRODUCTION 


Because of the current interest in thermal radiation effects during hyper- 
sonic flight, an investigation of the thermal radiation from ablation products in 
the stagnation region of blunt bodies is being conducted in ballistic ranges at 
the Ames Research Center. The present paper is a preliminary report of results 
obtained to date for three ablation materials in the velocity range of 4.5 to 
7-5 km/sec. The test technique consists of launching small diameter plastic or 
metal models at high velocities and measuring the resultant radiation with photo- 
electric devices. Previous publications describing this work are references 1, 

2, and 3* References 1 and 2 describe air radiation from the shock layer and 


reference 3 describes radiation from the wake. During the progress of these 
tests, it was noted that the spectra observed for a given test condition were a 
function of the model material, and, because of the ablation occurring on the 
front face of the model, these anomalies were attributed to foreign material 
radiating in the boundary layer. In certain instances, the intensity of the 
radiation presumed to be due to the foreign material was equal to or greater 
than that due to air. The question was raised as to the importance of the abla- 
tion product radiation to the total surface heating for full-scale vehicles. A 
study was therefore undertaken to quantitatively evaluate ablation product radia- 
tion as a source of radiative heating to vehicles in hypersonic flight. 

Prior work, both analytical and experimental, related to the radiation 
associated with the ablation process is described in references 4 and 5- The 
experiments consist of measurements of the radiation resulting from the ablation 
of plastic samples in a supersonic air stream heated by an arc discharge. The 
enthalpy range of the tests was equivalent to velocities from 1.5 to 5 km/sec. 

In the analysis of these data, the radiation was assumed to come from the solid 
surface of the material under test. 

The present investigation included both experimental work to measure the 
radiation from the ablation products at model scale and theoretical work to 
estimate the radiation from this source at full scale. The experiments were made 
with three ablating plastics and with aluminum (which did not ablate under the 
test conditions) at velocities from 4.5 to 7 .6 km/sec and at free-stream air 
densities from 0.02 to 0.2 amagat. 


SYMBOLS 

C boundary-layer-edge velocity gradient, sec" 1 

E radiant output per unit mass, watts/gm 

k ratio of universal gas constant to molecular weight 

K radiation constant, watts (cm) 3 ^ (gm) (°k) 

m mass injection rate per unit surface area, gm/cm 2 -sec 

M.W. molecular weight 

p pressure, gm/cm 2 

Pr Prandtl number 

q surface heating rate, watts/cm 2 

r radius from line of symmetry to the location s on the surface of 

vehicle, cm 

R body nose radius of vehicle, m 
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distance along "body surface, measured from stagnation point, cm 
Schmidt number 
temperature, °K 

velocity in s direction, km/sec 
velocity, km/sec 

mass fraction of ablation products in the boundary layer 
distance measured normal to surface, cm 
value of y where W = 0, cm 
transformed normal distance, p w dz = p dy 
boundary-layer thickness 

temperature difference across boundary layer, T e - T w , °K 
wavelength of peak response of the photodetector, microns 
universal gas constant 

exponential dependence of E on density 
exponential dependence of E on temperature 
density, gm/cm 3 

angle subtended by radius from axis to location s on surface 
value at the corner of the face (taken to be the sonic point) 

Subscripts 

ablation 

convection 

edge of boundary layer 

temperature 

wall 

ablation product mass fraction 


free stream 
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ablation products 


2 air 

o standard atmospheric conditions 


TESTS AND APPARATUS 


Small models were launched from a light -gas gun through air at selected 
pressures and the radiation from the model flow field was measured by photometric 
devices in the test section. The model launcher used is a two-stage, shock- 
heated, light-gas gun with a bore diameter of approximately 7 mm (refs. 1, 2, 3)* 
The test section is located 6 meters beyond the muzzle and is isolated from the 
muzzle blast by means of a blast deflector plate and a vented section. 

The main instrumentation (fig. 1) consisted of the spectral detectors for 
measuring radiation. In addition, two shadowgraph stations were equipped with 
electronic timing devices (for velocity and attitude measurements) and a Kerr 
cell shutter camera (for single frame photographs of luminous material near the 
model). The spectral detectors used in these experiments, one of which is shown 
schematically in figure 2, were multiplier phototubes in combination with narrow 
pass optical filters. The response of each combination was therefore limited to 
a narrow spectral range. By the use of eight such units, low resolution spectra 
of the radiation generated by the models were obtained. Shown in the table below 
are the characteristics of the phototube -filter combinations. The band passes 
listed are the widths between the half -maximum sensitivity points. 


Maximum sensitivity, 
wavelength 

Band pass 

Tube type 

Photoemissive 

surface 

0.235i_i 

0.049m 

RCA 

1P28 

S-5 

.265m 

.049m 

RCA 

1P28 

s-5 

• 312h 

.0561-1 

RCA 

1P28 

s-5 

• 369m 

.0381a 

Dumont 

6467 

S-ll 

. 460|_i 

.026p 

Dumont 

6467 

S-ll 

• 518m 

• 059m 

Dumont 

6467 

S-ll 

. 66lp 

•077m 

Dumont 

K1430 

S-l 

• 942p 

• 175m 

Dumont 

K1430 

S-l 


The slits, shown in figure 2, limit the field of view of the phototube to 
approximately 1 cm of flight path. The geometry of the slits permits the entire 
shock layer to be viewed, and the shock-layer and wake radiation to be dis- 
tinguished. 

Because the sensitivity of a multiplier phototube is strongly dependent on 
the high voltage supplied to the dynodes, the power supplies were carefully 
regulated and monitored. During many of the tests, neutral optical filters were 
used to avoid saturation and resulting nonlinear operation of the phototubes. 

All filters and phototubes were carefully calibrated by procedures described in 
appendix A. 
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Models of two different geometries were used, as shown in figure 3* In 
order to investigate the relationships between model material and ablation- 
product radiation, four materials were tested: aluminum, the polycarbonate 

lexan, polyethylene, and the General Electric Company Century Series plastic 124 
(hereinafter referred to as G.E. 124). Some materials which would not withstand 
launching stresses in the basic configuration were launched successfully when 
supported by a sabot. The sabots were made of lexan. For all the data in this 
report, the sabots were separated sufficiently from the model at the measuring 
station so that the model radiation could be monitored without interference from 
the sabot. 

The composition of the test air was monitored and no appreciable con- 
taminants were found. 


TEST RESULTS 


During the tests, the responses of the phototubes were displayed on 
oscilloscopes. These displays were photographically recorded and a complete set 
from a typical test is shown in figure 4. The wavelength of maximum sensitivity 
of each spectral detector is noted below each record. At wavelengths in the 
ultraviolet, the signal obtained consists of a single pulse which is interpreted 
as resulting from energy radiated from the shock layer. At the longer wave- 
lengths, the signal undergoes a distinct change in wave form; this change 
includes the development of a small rise immediately added to the initial pulse 
and the appearance of a second pulse due to energy radiated from the wake. The 
first rise of the initial pulse is from the nose region traversing the region of 
partial view (fig. 2). The further rise of this pulse occurs as the luminosity 
along the sides of the model comes into view. The relationship between phototube 
response and model luminosity is illustrated by the self -luminous photograph and 
the oscilloscope record in figure 5- The time scale in the oscilloscope record 
has been adjusted to match the distance scale of the photograph. The photograph 
was taken with the Kerr cell shutter camera at 50 nanosecond exposure time with 
film responsive in the wavelength region from 0.42 to 0.6[i. 

The basic data of the investigation which consist of measurements of the 
radiation from the nose region were obtained by reading the height of the first 
pulse and, by the method outlined in appendix A, converting to a radiant inten- 
sity of the source in watts/micron. The proper reading point for the height of 
the pulse is the end of the first rise, taken to be the entire nose region in 
full view, but before the radiation from the model sides is in full view. 

Typical results obtained by this procedure are shown in figure 6 for the 
four materials tested at three air densities. The short horizontal bars in the 
figure are the data reported in watts/micron radiated into 4rt steradians. The 
lengths of the bars are the band pass widths of the detectors as defined earlier. 
Lines have been faired through the data to approximate the spectra. 

The expected radiation from the air in the shock layer is shown superimposed 
on the data for the aluminum models. The prediction is taken from reference 6 
for air at equilibrium. The predicted level of radiation for air agrees well 
with the total shock-layer -radiation data for the aluminum models at the two 
higher air densities. The disagreement at the lowest density is attributed to 
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nonequilibrium air radiation (ref. 2). Some disagreement is noted, at all den- 
sities, however, at wavelengths greater than 0 . 5 m- • Whether this difference 
indicates a deficiency of the air -radiation prediction is not known at this time. 
Calculations made for the surface temperatures of the aluminum models indicated 
that melting temperature was never reached when the models passed the radiation 
detectors. Therefore, it is concluded that no aluminum ablation products were 
present and the spectra are representative of air. 

Based on the above, the spectra obtained with the aluminum models were 
subtracted from those obtained with the plastic models to remove the air- 
radiation contribution. (A slight adjustment was necessary to correct for the 
small difference in size between the aluminum models with sabots and the full- 
bore plastic models.) The results of this procedure are shown in figure 7 - 1 It 
is clear from this figure that a considerable quantity of radiation in excess of 
that expected from the air is observed with these materials, and the amount is 
strongly dependent on the material itself. 

The non-air radiation can be due to either the hot ablating surface or to 
the products of ablation in the boundary layer. If it is due to the former, then 
it is a cooling mechanism for the heat shield. If the latter is correct, then it 
constitutes a potential mechanism for additional heat input to the body. 

A simple test was devised to determine the origin of the non-air radiation. 
It consisted of launching models with effectively no projected nose area (i.e., 
right circular cylinders) and contrasting the measured radiation with that meas- 
ured from spherical faced models. If the non-air radiation were predominantly 
coming from the hot surface, the observed output would be severely attenuated 
because the projected area of the model surface viewed by the detectors would 
be essentially zero. Figure 8 gives results of these tests with lexan models 
together with results from one lexan model with the usual geometry chosen for 
similar test conditions. Flattening the model face increased the radiation 
observed in spite of the facts that less surface area was in view of the radia- 
tion detectors and the convective heating rate was lower. In addition, an effort 
was made to possibly attribute this radiation in part to a radiating surface in 
view of the detectors such as the corners which could ablate more rapidly than 
the stagnation area because of high convective heating. For the 6.72 km/sec test 
it was estimated that the corner rounded to give a projected area of not more 
than 10 “ 3 cm 2 during the flight time. Assuming a surface temperature equal to 
melting carbon, and an emissivity of 1 . 0 , a surface of this size could radiate 
only 5 percent of the observed amount. This shows that the non-air radiation 
emanated primarily from ablation products which are injected into the boundary 
layer. The same tests were conducted on G.E. 12k and polyethylene with analogous 
results . 


■'•Because of the small variation in velocity between various tests, a small 
velocity correction was made before subtracting the air radiation from the total 
spectrum obtained from the plastic materials. Reference 1 shows that over the 
range of these tests the air radiation is proportional to the eighth power of 
the velocity. This relationship was therefore used for the velocity correction. 
Generally, these corrections were less than 15 percent, but in a few instances 
were as high as 60 percent. 
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Another basic question was whether or not steady-state ablation was 
established for the plastic models in these tests. For the case of noncharring 
materials, the surface temperatures of the models were calculated as a function 
of flight-path distance for the test conditions of lowest heating rates. It was 
found that the model always traveled at least 20 times the distance required for 
the surface temperatures to reach that necessary for ablation, intimating that 
sufficient time is available to establish steady-state ablation. 

To calculate the quasi -steady-state conditions of the char -forming G.E. 124, 
a model was hypothesized as follows: Because of the very short time of flight 

(l millisecond) the char -spallation sequence was assumed to be in its initial 
cycle (ref. 7)* The char thickness was assumed to increase linearly along the 
flight path as the plastic pyrolized, and the preheat time was assumed short 
compared to the flight time. Char thicknesses were calculated and were typically 
of the order of 10“ 3 cm at the test section. Because the char was very thin, it 
was taken to have little effect on the ablation properties of the material, 
influencing only wall temperatures . 

The data of figure 7 were converted to average heating rates of the model 
face by integrating the spectra from 0.45b to l.lq and dividing half the integral 
by the surface area that the observed radiation falls on, namely, the curved sur- 
face of the model face. Because of the data scatter, intensity measurements 
were limited to spectral regions where the ablation product radiation was at 
least greater than the uncertainty of the measurements of the air radiation; 
hence, the lower limit of integration is 0.45 micron. It is to be noted, how- 
ever, that in some specific tests ablation product radiation below this limit 
appears to be measured. It is concluded that this additional apparent radiation 
is due to scatter in the data, and is not interpreted as ablation product radia- 
tion. The upper limit represents the long wavelength cut-off of the instrumenta- 
tion. Because of this restricted wavelength region, the resulting heating rates 
represent lower limits of the actual ablation product radiation heating rates . 
Allowance was made for that portion of the radiation hidden from the detector by 
the convex model face. Results of the integration are given in table I along 
with other parameters of interest for each test. It is apparent that under the 
test conditions, the amount of radiation from lexan ablation products was con- 
siderably more than radiation from ablation products of the other two materials. 


ANALYSIS 


It is not possible at the present state of knowledge to treat in complete 
detail the physical processes involved in ablation product radiation which would 
require taking into account specific chemical reactions, species present, etc. 
Instead, in the analysis which follows, consideration is given to gross boundary- 
layer processes which affect this radiation. The analysis relies on the present 
experiments and has as its objective the rational extrapolation of these experi- 
mental results to other air densities and larger vehicles. 
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Mathematical Model 


As was demonstrated in the last section, the non-air radiation measured 
during the present study could not have originated from the model surface; hence, 
the following analysis is based on the initial assumption that the radiant 
sources are ablated materials which are diffusing into the boundary layer over 
the model face. It is further assumed for simplicity that all important chemical 
reactions occur near the wall; subsequent reactions do not seriously alter the 
radiant power output. 

In general, one would expect the total power output per unit mass of foreign 
material (defined as unreacted ablation products and ablation products which have 
reacted with air) in the boundary layer to take the form E T = K T cp(T) where cp(T) 
is some pure function of the temperature. It is more convenient, however, to 
correlate the radiation on the basis of the mass loss from the model surface 
instead of the mass of foreign products in the boundary layer since it is very 
difficult to predict the extent of chemical reactions at the wall and, hence, the 
actual mass of foreign materials present in the boundary layer. However, these 
reactions may be sensitive to changes in density; therefore any correlation based 
on mass loss must provide for a density dependence. 

The following expression is then postulated as the power output per unit 
mass of ablation products (ablation product being defined as the material in the 
boundary layer which was originally part of the surface): 

E = KpM (l) 


where the effect of any chemical changes will appear especially in the magnitude 
of Note that if £ is zero, the density term in the above equation is 

removed and this infers that there are no chemical reactions present, or, if 
there are, they produce no change in the radiant energy output, or the reactions 
are not density dependent. Chemical reactions may also contribute to the magni- 
tude of however, no particular interpretation may be assigned to any specific 
value of | . 

The local heating, q^, due to radiation from a boundary layer, which is thin 
compared to its lateral dimensions, is obtained by taking half of the integral of 
the radiant output as a function of normal distance to the wall, namely. 


<3«. = y 


WEp dy 


( 2 ) 


where W is the mass fraction of ablation products in the boundary- layer mix- 
ture, p is the total gas density (air and ablation products), both being func- 
tions of y, and Y is the value of y when W = 0. 
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We introduce a transformed distance z 


which is defined by the equation 


p dz = p dy (3) 

W 

where p is the value of p at y = 0. 
w 

We now assume a linear dependence of W and T on this transformed variable: 


T e - T z 

— = 1 - £ 

T e - T w & T 


(4) 


_W _z_ 


(5) 


As used here, the temperature boundary -layer thickness, &j, is defined as that 
value of z where T/T e =0.99 and similarly the concentration boundary-layer 
thickness, 8^, is that value of z where W/W w = 0.01. 

These assumptions appear to be reasonable in the light of results obtained 
by Howe (ref. 8) and Roberts (ref. 9)- Using a detailed method to calculate the 
mass distribution of foreign species injected into the boundary layer, Howe found 
a nearly linear distribution of W with respect to r), which is a transformed 
distance related to z by a constant determined by vehicle geometry. Roberts, 
using the same linear assumptions as above, developed equations which describe 
the heat -transfer characteristics of a boundary layer with mass injection. He 
found that these equations gave results that agreed very well with boundary- 
layer heat transfer calculated by several exact methods . 

If one assumes the mixture to be a perfect gas, in order to use simple 
relationships between the state variables, the equation of state can be written 


p = (k lPl + k 2 p 2 )T 

Here k x is the gas constant divided by the average molecular weight for the 
ablation products and p is the partial density of the ablation products. The 
terms k 2 and p £ have analogous meanings for air. 

Since p = Wp and p £ = (l - W)p, the equation of state can be written in 
the form 


p = (k 3 W + k 2 )pT 


where 


k 3 = k x - k 2 
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Thus the total density may he -written as a function of the concentration as 


P 

(k 3 W + k 2 )T 


( 6 ) 


By substituting equations (l), ( 3 ), (4), (5), and (6) into (2) and retaining 
W as the variable of integration, we get 



2W w 



|^T e + AT 



1 


_W V 

W v 5 T/ 



(k 2 + k 3 W) 




W dW 


(7) 
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At this point we introduce the con- 
cept of a mass balance to determine an 
expression for in terms of ratios 

of boundary -layer thicknesses. The 
sketch of the forward surface of a vehi- 
cle with spherical geometry will be 
helpful in visualizing the following 
development . 

The mass injected into the boundary 
layer from s=Otos=s(s is the 
distance along the body surface measured 
from the stagnation point) is the inte- 
gral over the surface J m da, where m 

is the mass injection rate per unit sur- 
face area. The ablation material con- 
vected through a surface normal to the 
body at s = s is atrf* Wup dy evalu- 


ated at s. Here r is the radius from 
the line of symmetry to the plane s, and u is the gas velocity parallel to the 
surface. Noting that r = R sin cr and recalling that p dy = p v dz we may 
equate the two above expressions 


'o 

surface 


m da = (^2:rtR sin a p v 


&w 


Wu dz 

'o 7 at s 


Using equation (5), we write 



surface 


2jtR sin 0 


P Ww 
w w 



&W 




u dz 


J at 


s 


(8) 
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As before, we linearize u on z; u/u e = z/& u . Here &u is defined as that 
value of z when u/u e = 0 . 99 * Equation (8) may he integrated and written as 


/ m da = 
surface 

If the analysis is restricted to relatively blunt configurations, the convective 
heat transfer is not expected to vary more than 30 percent from that at the stag- 
nation point (ref. 10 ). To within this accuracy then, the mass injection rate 
may be considered a constant over the surface. The left side of equation ( 9 ) may 
then be approximated by 



s 

m da = [2jtR 2 (l - cos a)m] at g (10) 

surface 

We note also that, for blunt configurations, the boundary -layer edge 
velocity, u e , is a nearly linear function of s (ref. 11), increasing from 0 at 
the stagnation point to sonic velocity at the corner. The edge velocity may then 
be written 

u e = Cs = CRcr (11) 

where C is the velocity gradient. Thus, by utilizing equations ( 9 ), ( 10 ), 
and (ll) we can write the continuity equation: 



m 


1 - cos a 
a sin. a 


= 6 p w W ^ 



&W 


The trigonometric portion of the above expression is a slowly changing function 
of a, having a value of 0.500 at a = 0, and a value of 0.527 at a = 45°; it 
therefore may be replaced by l /2 without introducing significant error. The 
continuity equation then becomes 




( 12 ) 


This expression for Sy may then be used in equation (7) to obtain the heating 
equation in complete form: 


<1 


a 


f s v ^ w T , 

3 Kp^m /&u\ P 

2 W w 2 C W JJ 0 


T e + AT 



(k 2 + k^W)^ 




W dW 
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This equation can he rewritten 


. _ 3 Kp^m /&u i\ 1 P W 

= 2 C W W w £ J 0 



( k 2 + k3¥ ) 


W dW 


(13) 


In reference 9, Roberts has shown that the ratio of the concentration boundary- 
layer thickness to the velocity boundary -layer thickness is Sy/^ = (Sc w ) -0-3 
and since the Schmidt number for the mixture at the wall is estimated to lie 
between 0.8 and 2, it will be assumed, for the sake of simplicity, that by = & u . 
The ratio of the velocity boundary -layer thickness to the temperature boundary- 
layer thickness is also nearly 1 since, again, in reference 9, S^/bij = (Pr w ) 0,3 
and the Prandtl number is estimated to lie close to 1. With these two relation- 
ships we can write equation ( 13 ) as 


q., 


■a 


3 Kmp^ 1 

2 C W w 2 




( k 2 + k 2 W ) 




W dW 


(U) 


Recalling that k 3 is defined by k 3 = ki - k 2 , and k x and k 2 are ratios 
of the universal gas constant to the molecular weights, we can write the denomi- 
nator in the integral as 


(M.W. ). 


(m.w. ) , (m.w. ) , 


w 


A value of 29 is appropriate for (M.W. ) 2 and from reference 12, 36 seems to be 
a representative value 2 of (M.W. ) 1 . These values make it apparent that., k 3 W is 
small compared to k 2 . In fact, if the denominator is set equal to k 2 ^ and if 


T- 


I+i 


w 


is very small compared to T< 


l+i 


the heating equation may be integrated 


in closed form to obtain 


• _ 3 Km 5 

'la 2 C P e 


£+2 


(at ) 2 ( t - i + 1)(| -5+2) 


(15) 


To ascertain the error introduced hy this simplification, equation (l4) was 
machine -programmed and evaluated for the above -stated values of ’k 1 and kg. 


2 In reference 5, the charring ablation process was simulated experimentally 
and under certain conditions a value of approximately 5 'was obtained for the 
molecular weight of the ablation products. This value used in evaluating q^ 
from equation (l4) resulted in predicted heating rates only 8 percent lower than 
those obtained using a molecular weight of 36 . 
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When typical values for the other parameters were used, the value of q a given 
hy equation ( 15 ) was found to he within 5 percent of the value of q a given hy 
equation (14). 

The effects of changes in vehicle size on q^ are directly apparent from 
equation (15) because only m and C are sensitive to vehicle size. If the mass 
injection rate, m, is controlled primarily by convective heating, m is propor- 
tional to R' 1/2 and, if the sonic point is fixed somewhere near the corner of 
the face, l/C is proportional to R so that equation (15) reduces to q^ 
proportional to R 1/2 . 


Application to Test Data 


In order to use the test data to evaluate the empirical constants, equa- 
tion (15) is written in logarithmic form. Since the data cover a relatively 
small temperature range, | cannot be deduced directly without large uncertainty. 
Hence, an appropriate value will be assumed and the following equation solved 
for £ and K: 


log 


2q a C(AT) ; 

3mT e l+2 


= log K + £ log p 


- log(| - £ + 1) - log(| -5+2) 


Further, if the first two terms of the rapidly convergent series expansion for 
log (a + x) are used to expand the last two terms in the above equation, and if £ 
is anticipated to be much smaller than 2(| + 1), the above equation is written as 


log 


-2q a C(AT) 2 (| + 1)(£ + 2)' 


3mT, 


|+£ 


= log K + £ 


log p e + 0.434 (jA- 


I + 2. 


( 16 ) 


Thus, it is apparent from the above equation that if X x is represented as a 
function of X 2 , where 


X x = log 


-2Cq a (AT) 2 (£ + l)(| + 2)‘ 


3mT, 


| +2 


and 


X 2 = log p e + 0.434 (jA- + _L_) 

the resulting function is in slope -intercept form, where £ is the slope and 
log K is the intercept. 
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At this point it is appropriate to discuss "briefly the methods used for 
calculating the values of the parameters appearing in X x and X 2 . Parameters not 
dependent upon the material hut, rather, on the state of the air in the shock 
layer are the temperature and density at the edge of the boundary layer, T e and 
p e , and the velocity gradient along the surface at the edge of the boundary 
layer, C. The first two of these values, T e and p e , are taken to be those behind 
a normal shock wave for high temperature air (ref. 13)* The velocity gradient, 

C, was calculated assuming a linear velocity history along the edge of the 
boundary layer from the stagnation point to the sonic line. 

The mass injection rate, m, and the temperature difference across the 
boundary layer, AT, are dependent upon the properties of the ablating material. 
Reference 14 was used to calculate m for polyethylene and lexan tests; informa- 
tion was not available for calculating m for G.E. 124, so it was assumed to be 
like lexan in this respect. To calculate the temperature difference across the 
boundary layer, a knowledge of the surface temperature of the ablating material 
is required. Reference 15 presents the surface temperature of ablating poly- 
ethylene. The surface temperature of lexan was unknown but it was assumed that 
lexan would behave similarly to polyethylene . It should be pointed out that the 
temperature difference across the boundary layer will be of the order of 6000° K 
so that an error of 200 ° K in estimating the surface temperature will change the 
temperature difference by only 3 percent. The method used in reference 7 for 
calculating surface temperature of charring ablators was employed to calculate 
the surface temperature of G.E. 124. 

Before values of and X 2 can be calculated and plotted, %, the remaining 
empirical constant, must be determined. The range of test velocities and conse- 
quent variation in T e was not considered large enough to evaluate £ with 
precision. Reference 1 6 presents theoretical and empirical bases for assigning 
a value of 4 to this parameter. Since an uncertainty in | will not critically 
affect extrapolations of ablation product radiation heating rates as long as 
predictions are not made to velocities greatly different from those of the 
tests, a value of 4 was chosen, (in the evaluation of % from the lexan data 
the values of K and shown below, gave 3-5 ±4*) 

Each parameter required for evaluating X x and X 2 in equation (l6) was 
calculated for the tests and is listed in table I. Test values of Xx and X 2 
are shown in figure 9* The empirical constants, K and £, are then found to be 
as follows 

K S 


Lexan 

7x10' 8 (average) 0 

(average) 

G.E. 124 

2x10' 8 

.12 

Polyethylene 

6x10 " 7 

.66 


In this figure, the large symbols used for the lexan data represent an arithmetic 
average of the data at each of the three test densities. The actual test values 
are indicated by the smaller symbols. It may be seen that the possible value 
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of £ for lexan could range from - 0.18 to +0.22, with the value of 0 "being 
indicated by the average points. The range of possible values of £ for lexan 
indicates that the accuracy in determining £ by the present method is approxi- 
mately ±0.2. It should, therefore, be noted that £ for G.E. 12b could also 
vanish; however, the determined value of t, for polyethylene is almost certainly 
greater than zero. This is significant inasmuch as it indicates that the density 
dependence of the chemistry of ablation products is different for polyethylene 
than for lexan or G.E. 12k because of the chemical dissimilarity of polyethylene 
and the other two plastics. 


Application to Large Size Vehicle 


To assess the importance of ablation -product radiative heating for large 
size vehicles, equation ( 15 ) was evaluated for a vehicle with a 5 -meter nose 
radius; the convective heating only was used to calculate the mass injection 
rate, ignoring radiative heating from air or ablation products. If these radia- 
tive sources become important, they can greatly increase the ablation rate, since 
there is negligible blocking of radiation by ablation products in the boundary 
layer. 3 The flight conditions used in these calculations were an altitude of 
70 km and a velocity of 7 km/sec. The experimentally determined values of K 
and £ for each of the ablation materials tested were used in these calculations. 
The assumption is made in the analysis that t, is a constant over the extrapo- 
lated range, depending only on material. The actual full-scale heating may 
depart considerably from the derived model if, for example, radical changes in 
the chemistry of ablation product radiation occur in the extrapolated range. The 
calculated ablation product radiation is listed below together with the heating 
from air radiation and the cold wall convective heating. It is to be noted that 
for some materials the ablation product radiation can exceed the air radiation, 
and the resultant heating can be a substantial fraction of the cold wall convec- 
tive heating. Also to be noted is the wide range of heating rates, depending on 
ablator material. The ablation product radiation stated for lexan under these 
conditions is high enough to be self -aggravating as defined in the footnote in 
this section, if the surface absorptivity is greater than a few tenths. 

The local exponential dependence of the ablation product radiation on free- 
stream velocity and density was calculated and is included. 


3 It can be shown by writing equation ( 15 ) in difference equation form, and 
including the ablation product radiation to calculate the mass loss, that unless 

m) < h 

where a is the absorptivity of the surface, h is the heat of ablation, and 
(4a/" 1 ) is taken directly from equation ( 15 ), the process will he self -aggravating 
and the heating will approach limits imposed by available energy in the airstream. 
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4 , 

Material or heat source watts/cm 2 


Lexan 24 

G.E. 124 1.5 

Po ly ethylene . 0 43 

Air radiation .67 

C. W. convection 35 


Exponential dependence on 


Velocity 

2.6 

2.6 

3.0 


Density 

O .67 
• 77 
1.3 


The exponential dependence on velocity of ablation product radiation varies 
with velocity because of the varying relationship that exists between velocity 
and temperature of shock heated air. The exponential dependence on density is, 
as expected, strongly influenced by the value of £ for a specific material, 
since the boundary-layer-edge density is nearly a linear function of the free- 
stream density. 


CONCLUDING REMARKS 


In the foregoing analysis a source of vehicle heating not ordinarily 
considered was investigated - the thermal radiation from the gaseous products 
injected into the boundary layer by the ablation process. Tests conducted on 
small size plastic and metal models showed that important amounts of radiation 
can be contributed by these products. The experimental data were used to evalu- 
ate empirical constants in a derived ablation radiative heating equation. From 
the resulting equation predictions for heating from this source were made for a 
large vehicle. 

The results indicate that the ablation product radiative heating on a large 
vehicle can be a substantial fraction of the total heating, depending critically 
on the flight conditions and heat shield material; in some instances, this heat- 
ing exceeds the air radiative heating. 

The total heating is not necessarily the sum of all heating components 
presented. Each process, lowering the available energy in the shock layer, tends 
to moderate the other; hence, the total heating should be expected to be some- 
what less than the component sum. 

Since this study was a preliminary investigation of the ablation product 
radiative heating, and the extrapolations on size and density were made to 
several orders of magnitude beyond the test data, further investigations are 
necessary, particularly at lower densities, at higher velocities, and with larger 
size models. Materials representative of current ablative heat -shield material 
should be tested since it has been shown that the ablation product radiative 
heating is very strongly dependent on material. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., May 28, 1963 
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APPENDIX A 


CALIBRATION OF SPECTRAL RADIATION DETECTORS 
ABSOLUTE SPECTRAL CALIBRATION 


The absolute spectral response of a detector is determined by the narrow 
pass optical filter, the spectral response of the photocathode, the amplification 
of the multiplier phototube, and the construction of the voltage divider circuit. 

Relative spectral response, R(A) , of a spectral detector was defined by the 
formula 


= R'(a)t(a) 

[R' (A)T(A) ]max 


where R' (A) is the relative spectral response of the unfiltered photocathode and 
T(A) is the spectral transmission of the narrow pass filter. 


The absolute sensitivity of a detector is obtained by irradiating the 
detector from a unit distance with a radiant source, W (A), of known absolute 
spectral intensity. The voltage output is then 


e = K 



W(A)R(A)dA 


where K is the absolute sensitivity of the detector to radiant energy and has 
units of watts/volt. Since e, W(A), and R(A) are known, we can calculate K. 
When recording the radiant output from an unknown source, J'(A), we make the 
following approximation 



R(A)dA 


J'(A)R(A)dA 


where A 0 is the center of the interval where R(A) is nonzero. An average 
radiant flux for this interval is then represented by j(A Q ). 

The subsequent output voltage is then 


e = K 



J'(A)R(A)dA 


= KJ(A 0 ) 



R(A)dA 


and J(A 0 ) follows directly. 


IT 


A more general treatment of calibration of photodetectors will he found in 
reference 17 . 


PHOTOCATHODE SATURATION 


Each detector was calibrated to determine the maximum output voltage it 
will produce before the output becomes a nonlinear function of the incident 
radiation. 

A short -duration, high- intensity light pulse which simulated in duration 
and shape the radiation measured during the tests was imaged through a lens on 
the photocathode. An iris in the lens varied the intensity. This configuration 
was then used to measure the apparent transmission of a neutral pass optical 
filter of 10- to 20 -percent transmission. The apparent transmission is given by 

Output voltage with filter in light path ^in 

Apparent transmission = — — — ^ ~r = U 

Output voltage without filter in light path v 0 \xt 

When the light intensity is raised to the extent that the detector is saturated, 
V ou t is no longer a linear function of the incident radiation. Since V ou t is 
from 5 to 10 times greater than V^n, it was possible to cause V ou t to be a 
measurable amount nonlinear before Vin. was affected. This saturation was then 
noted as an abrupt change in the apparent transmission and V ou -t was noted as 
the limiting output voltage. 

The above was repeated for various values of the high voltage supplied to 
the divider circuit and the corresponding limiting voltages were noted. 

Neutral pass optical filters could then be used during the tests to assure 
that the detectors were used in this linear range. 
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TABLE I.- SUMMARY OF TEST PARAMETERS AND RESULTS 


Material 

V*, 

km/sec 

Poo/p o 

9a> 0 

watts / cm 2 

A 

gm/cm 2 sec 

le , 

°K 

°K 

4 C , 

watts /cm 2 

Aluminum 

6.47 

0.02025 











6.22 

.0800 

— 

— 

— 

— 

— 


5-93 

.1974 

— 

— 

— 

— 

— 

Lexan 

6. 56 

.02025 

21.5 

0.887 

7, 450 

1,000 

11, 600 


6.82 

.02025 

36 

.9435 

7,660 

1,000 

13,200 


5.90 

. 02025 

14.5 -17.5 

.78 

6,950 

1,000 

8,400 


7.62 

. 02025 

32.4 

1.08 

8,290 

1,000 

18,100 


6.08 

.0800 

69 

1.587 

7, 570 

1,000 

18, 400 


5.75 

.0800 

37 

1.463 

7, 228 

1,000 

15, 400 


4.54 

.1974 

25.2 

1.556 

5,700 

1,000 

12,000 


5.88 

.1974 

83 

2.386 

7,710 

1,000 

26, 200 


6.72 

.1974 

190 

2.873 

8, 594 

1,000 

39,000 

G.E. 124 

6.83 

.02025 

5.64 

.934 

7,665 

1,550 

13,000 


6.125 

.0800 

17.0 

1.60 

7, 560 

1,890 



6.275 

.1974 

29 - 32.6 

2.61 

8,120 

2,220 


Polyethylene 

6.89 

. 02025 

2.1 

.816 

7,710 

1,000 



6.08 

.0800 

8.3 

1.326 

7,560 

1,000 



5.88 

.1974 

21.2 

1.961 

7,680 

1,000 



21 











22 



23 


Figure 1.- Schematic representation of test section. 
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~K Q = 0 . 460|a A 0 = 0.5l8(a A 0 = 0.66ln A 0 = 0.942ia 

Figure 4 .- Set of phototube signals recorded during a typical test. Sweep rate is 5 microseconds 

per large division. 


ro 
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Figure 5-- Self-luminous photograph and oscilloscope record showing relationship 
between luminosity about model and phototube signal. 
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Radiation, watts/micron 


Polyethylene 
V a0 = 5.88 km/sec 


G.E. 124 

V <0 = 6.275 km/sec 


4 .6 .8 


.6 .8 1.0 
Wavelength, microns 


(a) <J P 0 * °' 20 



Lexan 

Voo= 5.88 km /sec 


Polyethylene 

=6.08 km /sec 


G.E. 124 
V =6.125 km/sec 


' Lexan 

V w =6.08 km /sec 

i i i i i l i I 


.6 .8 1.0 1.2 A 

\ 

.6 .8 1.0 1.2 

4 .6 .8 



Wavelength, microns 



Cb) 

p„/p 0 =0-08 


Polyethylene 
V = 6.89 km /sec 

1 - 

G.E. 124 
V = 6.83 km/sec 

““ Lexan 

V =6.56 km /sec 


Wavelength, microns 

M p > 0 - °- 02 

Figure 7.- Spectra of ablation product radiation from the shock layer of models 





Total spectral energy, watts /micron 



Figure 8.- Comparison of spectra from right circular cylindrical and basic 

model configurations. 






